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ABSTRACT 

We use extensive measurements of the cluster A1689 to assess the expected similarity in the dynamics 
of galaxies and dark matter (DM) in their motion as collisionless 'particles' in the cluster gravitational 
potential. To do so we derive the radial profile of the specific kinetic energy of the cluster galaxies from 
the Jeans equation and observational data. Assuming that the specific kinetic energies of galaxies and 
DM are roughly equal, we obtain the mean value of the DM velocity anisotropy parameter, and the 
DM density profile. Since this deduced profile has a scale radius that is higher than inferred from 
lensing observations, we tested the validity of the assumption by repeating the analysis using results of 
simulations for the profile of the DM velocity anisotropy. Results of both analyses indicate a significant 
difference between the kinematics of galaxies and DM within r < 0.3r v i r . This finding is reflected 
also in the shape of the galaxy number density profile, which flattens markedly with respect to the 
steadily rising DM profile at small radii. Thus, r ~ 0.3r v i r seems to be a transition region interior to 
which collisional effects significantly modify the dynamical properties of the galaxy population with 
respect to those of DM in A1689 

Subject headings: clusters: galaxies - clusters: individual: A1689 



1. INTRODUCTION 

The properties of dark matter (DM), which is be- 
lieved to be mostly cold and collisionless, have been ex- 
tensively explored in dynamical simulations. A much 
touted result from these simulations is the "universal" 
density profile (Navarro, Frenk, & White 1997, hereafter 
NFW; Moore et al. 1998) of galaxies and galaxy clusters. 
As the main mass constituent of galaxy clusters, DM 
largely self-gravitates and dominates the hydrodynam- 
ics of intracluster (IC) gas and the dynamics of member 
galaxies. Cluster DM density profiles deduced from X- 
ray observations (Pointecouteau, Arnaud, & Pratt 2005; 
Vikhlinin ct al. 2006; Schmidt & Allen 2007; Arnaud, 
Pointecouteau, & Pratt 2008), galaxy velocity distribu- 
tion (Diaferio, Geller, & Rines 2005), SZ measurements 
(Atrio-Barandela et al. 2008), and strong and weak lens- 
ing measurements (Broadhurst et al. 2005a, hereafter 
B05a; Broadhurst et al. 2005b, hereafter B05b; Limousin 
et al. 2007; Medezinski et al. 2007; Lemze et al. 2008, 
hereafter L08; Broadhurst et al. 2008; Zitrin et al. 2009, 
2010; Umetsu et al. 2010) are broadly claimed to be con- 
sistent with NFW profiles. However, the shape of the 
profile in the inner cluster region where it is predicted to 
have a characteristic radial slope of —1 has been deduced 
to be shallower in some studies (Kravtsov et al. 1998; 
Ettori et al. 2002; Sanderson et al. 2004) and steeper, 
around -1.5, in others (Fukushige & Makino 1997, 2001, 
2003; Moore et al. 1999; Ghigna et al. 2000; Klypin et 
al. 2001; Navarro et al. 2004; Limousin et al. 2008). 
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A more complete description of the clustering proper- 
ties of DM necessitates characterization of its phase space 
distribution. Whereas the DM density profile can be de- 
termined directly from analysis of X-ray, lensing, and SZ 
measurements, deducing the DM velocity distribution is 
considerably more challenging. On the other hand, the 
dynamical properties of member galaxies can be studied 
directly in terms of their density and velocity profiles, 
including a study of the radial behavior of the velocity 
dispersion. Also, the location of the velocity 'caustics' 
can now be studied both in individual massive clusters 
where the data quality is high (Lemze et al. 2009, here- 
after L09) and in composite surveys for which only lower 
quality measurements are available but statistical results 
can be derived (e.g., Biviano & Girardi 2003). 

The dynamical evolution of the galaxy population in 
a cluster is presumed to be largely collisionless follow- 
ing an initial phase of mean-field ('violent') relaxation of 
the main sub-cluster progenitors that merged to form the 
cluster. This collisionless behavior is expected particu- 
larly outside the central cluster region. As such, the basic 
dynamical characteristics of cluster galaxies are expected 
to resemble those of DM, which is strictly collisionless. 
For example, in the cluster 1E0657-56 - the "bullet clus- 
ter" - despite a recent collision of two massive clusters, 
the spatial distributions of DM and galaxies are quite 
similar. A conical shaped shock front is visible indicat- 
ing two clusters have passed through each other with 
an obvious collisionally merged gas distribution, but the 
galaxies and the lensing mass are largely intact, implying 
straightforwardly that the DM and galaxies are collision- 
less (Markevitch et al. 2002; Clowe et al. 2004; Bradac 
et al.2006). In the colliding cluster A520, on the other 
hand, a massive dark core is claimed to coincide with 
the central X-ray emission peak, but the region is largely 
devoid of galaxies (Mahdavi et al. 2008), though this de- 
pends on the way the weak lensing analysis is formulated 
(Okabe & Umetsu 2008). 
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With increasingly extensive and precise data, such as 
we have acquired for A1689, it is now possible to assess 
the collisionless nature of galaxies and DM by measuring 
the degree of consistency between the measured galaxy 
and DM density profiles and the profile of the DM veloc- 
ity anisotropy. Af 689 seems well relaxed (with possibly a 
small deviation from a relaxed state; Andersson & Made- 
jski 2004). It has a centrally located cD galaxy, and an 
X-ray emission region that is spherically symmetric (Xu 
& Wu 2002; L08; Riemer-Sorensen et al. 2008). The 
cluster has well-defined galaxy velocity caustics with no 
major infall of matter close to the virial radius (L09), 
and only a low level of substructure (Broadhurst et al. 
2005a,b; Umetsu & Broadhurst 2008). 

We have previously determined the DM and gas den- 
sity profiles in AI689 from a combined analysis of lens- 
ing and X-ray measurements (L08), using an approach 
that we refer to as model- independent, since we did not 
assume particular functional forms for the profiles. Ad- 
ditionally, we extended the analysis by including photo- 
metric and spectroscopic measurements of a very large 
number of galaxies in the AI689 field, from which we de- 
duced the positions and radial velocities of 4,76^11 cluster 
members. These results made it possible to deduce the 
galaxy velocity anisotropy profile, which was found to 
exhibit the expected behavior, varying between predom- 
inantly radial orbits at large radii towards more tangen- 
tial orbits near the center (L09). 

Here we show that with the above information we can 
infer the DM velocity anisotropy which we compare with 
the galaxy velocity anisotropy profile and with results 
from simulations. The paper is organized as follows; in 
§ [2] we describe our method for determining the DM ve- 
locity anisotropy. In §[3] we compare the DM and galaxy 
density profiles (§ !3.2p . derive the DM velocity anisotropy 
and compare it with the galaxy velocity anisotropy pro- 
file and results from simulations (§ 13. ip . and estimate 
the collisionless profile of cluster galaxies (§ 13. 3|) . We 
conclude with a summary and discussion in § [4j 

2. METHODOLOGY 

In this section we present the procedure for deriving 
the DM velocity anisotropy using results from our pre- 
vious analyses of the galaxy dynamics (L09) and the 
total mass density profile (L08) of A1689. In L08 we 
combined lensing and X-ray measurements to determine 
model-independent profiles of the gas and total mass den- 
sity profiles (i.e., without assuming particular functional 
forms) . In the second stage of the work (L09) the galaxy 
surface number density and the projected velocity disper- 
sion were included and analyzed using the Jeans equa- 
tion, from which we obtained profiles of the 3D galaxy 
number density and the galaxy velocity anisotropy. 

The dynamics of a collisionless gas are governed by the 
Jeans equation (Binney & Tremaine 1987) 

Id, 2 , 2&a? r GM 

^0^) + ^— = — , a) 

where i =DM, gal, and p-i(r) is the density of element i. 
The velocity anisotropy profile /3j(r) is 



where a r (r) is the radial velocity dispersion, and at(r) = 
C0 ( r ) = °^( r ) is t ne (ID) transverse velocity dispersion. 
Using eq.[l]for the galaxies, the degeneracy between <7 ga i ir 
and /3 ga i can be removed with sufficient spectroscopic 
data (L09). This procedure is obviously irrelevant in the 
case of DM, for which we have to adopt an alternative 
approach. 

The orbit of a test particle in a collisionless gravita- 
tional system is independent of the particle mass. This 
would presumably imply that once hydrostatic equilib- 
rium is attained, most likely as a result mixing and mean 
field relaxation, DM and galaxies should have the same 
mean specific kinetic energy, i.e., 

°DM,tot( r ) = °gal,tot( r ) i (3) 

where 

^totW = 4r( r )+4e( r )+44 r ) = <r(r) (3 - 2A(r)) . 

(4) 

Additionally, it is expected that the total specific kinetic 
energy of DM particles is proportional to that of the gas 
(e.g., Mahdavi 2001; Host et al. 2008), T oc cr t 2 ot . Ob- 
servational evidence for this scaling relation comes from 
combined X-ray and optical observations of groups and 
clusters for which the mean emission-weighted gas tem- 
perature scales roughly as the second power of the total 
galaxy velocity dispersion (Mulchaey & Zabludoff 1998; 
Xue & Wu 2000). 

The temperature and density profiles of IC gas can 
be deduced from X-ray spectral and surface brightness 
measurements. These profiles can then be used to de- 
termine the total mass distribution from a solution of 
the hydrostatic equilibrium equation. We have recently 
developed a model-independent joint lensing/X-ray anal- 
ysis procedure (L08) to examine the consistency of X-ray 
temperature and emission profiles with the lensing based 
mass profile, finding that the cluster temperature profile 
is systematically ~ 30-40% lower than expected when 
solving the equation of hydrostatic equilibrium using the 
lensing-deduced mass profile. This discrepancy may re- 
flect in part the ambiguity in deriving 3D temperatures 
from projected spectral X-ray data, stemming from the 
sampling of a range of gas temperatures along any given 
line of sight (Mazzotta et al. 2004; Vikhlinin 2006). This 
could also be partly related to the small-scale structure of 
the gas [possibly including relatively dense cooler clouds 
found in simulations (Kawahara et al. 2007)] which may 
result in a significant downward bias in temperature es- 
timates from spectral X-ray observations. The inferred 
temperature is also sensitive to instrumental effects, as 
has recently been deduced in the analysis of Chandra ob- 
servations of A1689 (Peng et al. 2009). Other reasons for 
temperature biases can be deviations from equilibrium 
and non-thermal pressure (Molnar et al. 2010). Even in 
the best-case scenario the temperature can be used as a 
reliable tracer for t7 t ot only in the inner part of the clus- 
ter, where hydrostatic equilibrium applies, and also only 
at radii larger than about 0.1r v i r , since at smaller radii 
the specific energy of the gas and DM may be different 
(Rasia, Tormen, & Moscardini 2004). 

For our analysis we assume an NFW-like profile for the 
DM density 



(2) 



PDMMocKr/rsXl + r/OT 1 , (5) 
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and take the DM velocity anisotropy to be constant, i.e., 
/3dm = const, since the data is not sufficient to mean- 
ingfully constrain more than one free parameter in this 
quantity. We then relate our model parameters to the 
DM radial velocity dispersion, <Tx>M,r> using the Jeans 
equation, where the (total) mass profile is taken from 
L08. Here we allow for the possibility of a difference be- 
tween the total density profile (directly measured by lens- 
ing) and the profile of just the DM. From the DM velocity 
anisotropy and radial velocity dispersion we can then de- 
duce the DM total specific kinetic energy. Best-fit values 
of Pum(t) and /3dm are determined by fitting the DM 
total specific kinetic energy to the galaxy total specific 
kinetic energy. More specifically, /3dm and pdm(0 are 
determined by minimizing x 2 = XL V T ( r i) ' C ' ^( r «)) 
where V(r<) = cr DMtot (ri) - ^gai,tot( r <) and C is the co- 
variance matrix of the measured <7g al tot (^i)- The galaxy 
total specific kinetic energy itself is not a direct measure- 
ment. It was derived as shown in eq. 2] from fgai.r and 
Pgsdir). Since the values of er g ai,tot at the various radial 
positions are derived from underlying parameterized 
expressions (as given below for /3 ga i(r)), their uncertain- 
ties are correlated. 

The derivation of /3 ga i(r) and cr ga i,tot(?"i) is carried out 
by following the same procedure as in L09, except for 
allowing a higher freedom of /3 ga i(r) at large radii. N- 
body simulations for a variety of cosmologies show that 
the velocity anisotropy has a nearly universal radial pro- 
file (Cole & Lacey 1996; Carlberg et al. 1997). In accord 
with the work of L09, /3 ga i(r) is taken to have the follow- 
ing form: 

/3 gal (r) = (Aj + Poo) (r ^)2 + I - A> , (6) 

where we note that /3oo = 1 was adopted by L09. Now, 
the total number of cr ga i ;tot (r 2 ;) bins cannot exceed the 
total number of free parameters in the expressions for 
<7gal,r and /3 ga i(r), which is 6, since a larger number would 
cause a complete degeneracy among the various values of 
c g ai,tot(^i)- Even taking 6 bins resulted in an unphysi- 
cal correlation matrix (i.e., one having a negative eigen- 
value), which still indicates a near-degeneracy. There- 
fore we adopted 5 radial bins of cr ga i,tot(?'i)) which was 
the maximum number for which degeneracy is not sig- 
nificant and error estimates are reasonable. 

3. RESULTS 
3.1. Velocity anisotropy profiles 

The DM velocity anisotropy, /3dm j was determined as 
described in § [2J with the constraint that the DM total 
specific kinetic energy must satisfy er DM tot > 0. The re- 
sulting acceptable fit, with x 2 /dof = 3.5/(5 — 3), is shown 
in figure [TJ The best-fit parameters of the analytical 
expressions for the DM density and velocity anisotropy 
profiles are given in table [1] It is important to note that 
while the errors on the parameters are rather large, in 
this fit we have constrained the DM parameters based 
only on the fit to equation [31 without assuming the DM 
density profile to be similar to the total mass density 
profile. Thus, the results allow a largely independent 
comparison between the consequences of assuming equa- 
tion [3] and the results of other observational probes of 



TABLE 1 

The values of the parameters of DM density and velocity 

ANISOTROPY. THE ERRORS ARE 1-cr CONFIDENCE. 
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r a [h _1 kpc] 
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Fig. 1. — Profile of the total specific kinetic energy of the galaxies 
(blue circles and 1— a errorbars) and the fitted total specific kinetic 
energy of the DM (red line). 

the cluster. 

The corresponding galaxy and DM velocity anisotropy 
profiles are plotted in figure [2] (top panel) together with 
their respective 1— a uncertainties. In figure [2] (lower 
panel) we compare the derived DM velocity anisotropy 
value to the profile derived from simulations. The current 
data only allow us to determine an overall, typical value 
of /3dm m the cluster, and given the rather large uncer- 
tainties, there is fair agreement between this value and 
the typical value of /3 ga i in this cluster, and also with the 
typical values of /3dm seen in simulated clusters. Note 
from the figure that while we allowed /3oo to be a free 
parameter for the galaxies, the best-fit value came out 
quite close to unity. 

We briefly describe the derivation of the DM veloc- 
ity anisotropy profile from results of an ENZO simula- 
tion (for a more complete description of this particular 
simulation, please refer to Hallman et al. 2007). Since 
A1689 is a moderately-distant (z=0.183), high mass clus- 
ter, M vir - 1.5 ■ 10 15 h" 1 M (Broadhurst et al. 2005a; 
Oguri et al. 2005; Limousin et al. 2007; L08; Umetsu 
& Broadhurst 2008; L09; Umetsu et al. 2009; Corless 
et al. 2009; Coe et al. 2010), we derive the DM veloc- 
ity anisotropy profiles for a sample of high-mass clus- 
ters with M vir > 10 15 h a7 M at z = 0.2. The clus- 
ters were drawn from a cosmological adaptive mesh re- 
finement (AMR) simulation performed with the ENZO 
code developed by Bryan & Norman (1997) and Norman 
& Bryan (1999) assuming a spatially flat CDM model 
(very similar to the concordance model). The AMR sim- 
ulation assumed adiabatic gas dynamics (i.e., neither ra- 
diative heating, cooling, star formation or feedback were 
included). The box size was 512 h _1 Mpc comoving on 
a side with 512 3 DM particles, and DM mass resolution 
of about 10 11 hp 7 M . The root grid contained 512 3 
grid cells, and the grid was refined by a factor of two, 
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Fig. 2. — Top panel: comparison between /9dm and /3 ga i. Shown 
are the velocity anisotropy profiles of galaxies (blue dotted lines) 
and DM (black solid lines). Bottom panel: comparison between the 
DM velocity anisotropy inferred from data (black solid lines) and 
that derived from simulations (blue dot-dashed lines). For each set 
the central line shows the best fit and the other two lines show the 
± 1— <j uncertainty range. 

up to seven levels, providing a spatial resolution of 6.5 
kpc h -1 at z = 0.2. For each halo out of the total num- 
ber of 51, we derive the DM velocity anisotropy profile, 
which is averaged over all halos, with an uncertainty that 
is taken to be the standard deviation as determined by 
Lemze et al. (2010, in preparation). The derived profile 
is in agreement with those found in other studies (e.g., 
Colin, Klypin, & Kravtsov 2000; DM04; Mamon & Lokas 
2005; Valdarnini 2006). 

3.2. Density profiles 

To further explore phase space occupation of DM and 
galaxies, we compare the deduced DM density profile to 
the galaxy and to the total density profile. The shape of 
the galaxy density profile is represented with a core, and 
that of the total density with either a model- independent 
profile, a 'universal' NFW form, or a cored form. In 
figure |3] we compare our results for these profiles, i.e., 
we compare the DM density from our current fitting to 
equation (green dot-dashed curves showing the best-fit 
result as well as the l—o~ uncertainty region) to our previ- 
ous results for the galaxy density (blue solid curves) and 
the total mass profile, whose various versions are shown 
by the points with error bars (model-independent fit), 
dashed black curve (NFW), and dotted red curve (core). 
In order to include all these in the same figure, and allow 
a comparison of the relative shapes, we arbitrarily scaled 
the profiles so that they match at 700 h _1 kpc ~ ir v ; r 




10 

Radius [h~ 1 kpc] 

Fig. 3. — The galaxy density (blue solid curve with upper and 
lower solid curves marking the 1— a uncertainty) is compared to the 
deduced DM density (green dot-dashed curves showing the best-fit 
and and 1-<t uncertainty) and to the total matter density for the 
model-independent fit (points with error bars), the NFW fit (black 
dashed curve) , or the core fit (red dotted curve) . We show relative 
density profiles, all scaled to match at 700 h _1 kpc (except for the 
galaxy 1— <r lines). The left and right black vertical lines indicate 
gfvit and r v i r , respectively. Note that the galaxy density line at 
low radii is an extrapolation due to lack of data in this region (and 
therefore no error bars are shown). 



(with arbitrary units in the j/-axis). 

While the DM scale radius derived here from the ve- 
locity dispersion fit has a large error, its best-fit value 
is significantly higher than that derived previously (L08, 
table 4) for the total (mostly DM) density. This is re- 
flected by the slower rise of the DM density at small radii 
compared to the total density (figure [3]). This may indi- 
cate that our assumption of eq. [3] invalid at small radii. 

Another perspective on the assumption that DM fol- 
lows the same dynamics as the galaxies can be seen from 
the steepness of the directly measured density profiles. 
To assess the steepness of the fitted profiles we plot in 
figure 2] the distribution of their power-law indices, 
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^logL°(r)] 
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(7) 



The power-law index of the galaxy profile is shown by 
the blue lines, and that of the total mass profile by the 
black (NFW) and red (core) lines. We did not plot also 
the power-law slope of the mo del- independent fit, since 
it is clear from figure |3] that it would look very similar to 
the results for the NFW and core profiles. Since the DM 
profile should be rather similar to the total mass profile, 
the conclusion from these figures here again is that the 
galaxies and the DM have consistent density profiles for 
r > ir v ir, but the profiles are significantly different at 
smaller radii. 

3.3. The collisionless profile of cluster galaxies 

In § 13.11 we derived the DM velocity anisotropy profile 
assuming that both galaxies and DM are collisionless, 
finding that the best-fit value of DM density scale radius 
is higher than inferred from direct lensing observations 
for the total mass (L08). As was mentioned above, this 
may indicate that eq. [3] is not valid at all radii. To quan- 
tify differences between <Tp M tot and <7g al tot we use the 
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Fig. 4. — Top panel: power-law indices of the galaxy and total 
mass profiles. The index of the galaxy profile is shown by the 
blue solid line, with upper and lower solid lines indicating 1-cr 
uncertainty, and that of the total mass is shown by the black dashed 
(NFW) and dotted red (core) lines, with 1— a uncertainties. Bottom 
panel: ratio between the power-law index of the galaxy density to 
that of the total matter (NFW - black dashed with upper and 1-cr 
uncertainty; core - red solid with 1-cr uncertainty). The left and 



right black vertical lines mark 



and r v i r , respectively. 
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'DM, tot 
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(r) 



(8) 



which we expect to be very close to unity if both compo- 
nents are fully collisionless (or if both deviate comparably 
from the purely collisionless limit). 

We have thus far assumed that f co \\(r) must equal 
unity and tried to obtain acceptable fits under this as- 
sumption, deriving other results in the process. Here 
we try a different approach, where we start by assum- 
ing that the DM velocity anisotropy profile in A1689 
matches the profile derived from simulations (see § 13. ip . 
We used our previously-determined DM mass density, 
Pdm — jOtot — /°gas, derived by assuming profiles for the 
total and gas density profiles. Two kinds of profiles were 
assumed, a model-independent and a model-dependent 
one, where in the model-dependent case we assumed an 
NFW and double beta model for the total and gas den- 
sity profiles, respectively. After assuming a profile for 
the total and gas density profile, we can then solve the 
Jeans equation to determine <7p M tot . In figures [5] and [6] 




500 1000 1500 

Radius [rf 1 kpc] 

Fig. 5. — Profiles of / co n (top panel, blue solid curve) and the 
velocity bias b (lower panel, blue solid curve) taking the total and 
gas density model independent profiles from LOS along with their 
1— a uncertainty regions (marked by blue dot-dashed and dashed 
lines, for the case when the uncertainty of DM velocity anisotropy 
from simulations is included, and the case when this uncertainty 
is not included, respectively). The vertical black line marks jr v i r , 
and the horizontal dotted line indicates the expected value if both 
DM and galaxies are purely collisionless. 

we show the resulting ratio / co ii(r), or equivalently the 
corresponding velocity bias of the galaxies relative to the 
DM 

b(r) = W/coiiW , (9) 

for the model-independent and model dependent pro- 
files, respectively. Our results are consistent with both 
the galaxies and the DM being purely collisionless at 
r > l^vir, but there is some evidence for a significant 
deviation from this limit at smaller radii (except for the 
very smallest radii, where the observational constraints 
are weak). 

4. DISCUSSION 

The work reported here is a continuation of our com- 
prehensive study of the dynamical properties of DM and 
galaxies, and the hydrodynamical properties of IC gas 
in the well-observed cluster A1689 (L08, L09). In L09 
we derived the galaxy density and velocity distributions, 
from which we deduced the specific kinetic energy of the 
galaxies. In the work reported here we assumed that if 
DM and galaxies are fully collisionless they should have 
the same average specific kinetic energy (as manifested 
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Fig. 6. — The same as figure \E\ except that here we take the total 
and gas density profiles to be NFW and double beta model (as 
described in L08), respectively. 

in eq. [3]) . Using the mass profile which was previously de- 
rived in L08, together with the Jeans equation, we fitted 
the DM to the galaxy specific kinetic energy and deter- 
mined the DM density and the DM velocity anisotropy. 

For the DM density we obtained a best-fit value for 
the power law index at large radii, a + 1 , which is higher 
(i.e., steeper) than the NFW value of 3 by about la. 
The best-fit value of the scale radius is ~ 1.9a higher 
than the scale radius of the total density This was the 
first of several indications that the dynamical state of 
the galaxies and the DM may differ at small radii, below 
~ ^r V i r . In particular, our assumption of eq. [3] resulted 
in a derived DM density profile at small radii that was 
inconsistent with (specifically, shallower than) that we 
have previously measured for the total matter density. 

In order to quantify the possible difference between 
°dm tot an d cr| a i tot * we use d another approach to deter- 
mine their ratio / co ii- Specifically, we adopted the profile 
of /3dm deduced from simulations and our derived DM 
mass density (from L08), together with the Jeans equa- 
tion. The deduced profile of / co n can be interpreted as a 
measure of how collisionless are the galaxies with respect 
to the DM (which is equivalently expressed also by b, the 
velocity bias of eq.|H]). Note that in some simulation work 
the velocity bias is defined in terms of the velocity disper- 
sions of DM subhalos (rather than galaxies) compared to 
the DM background. We find that / co n < 1 (at > la) at 



some radii below ^r v ; r , but that it is consistent with unity 
at larger radii. This implies that at r > ^ r vir we can in- 
deed assume the validity of eq. [3j but at r ~ 450 h _1 kpc 
we deduce b s» 1.5lto"i an d b « 1.3^;!, for the model- 
independent and model dependent profiles, respectively, 
an indication that b increases from large radii towards 
the center. This is in accord with results from simula- 
tions which found b to vary from slightly higher than 
one (by ~ 10%) at large radii (Ghigna et al. 2000; Colin, 
Klypin, & Kravtsov 2000; Diemand, Moore, & Stadel, 
2004, hereafter DMS04; Gill et al. 2004) to b ~ 1.3 in the 
central region (r < 0.3r vir ) (Ghigna et al. 2000; DMS04). 

We note that there is some uncertainty regarding the 
dependence of this bias on the features of the subhalos 
or galaxies considered; e.g., whereas DMS04 found the 
bias to be independent of the subhalo mass, Ghigna et 
al. (2000) claim that there is no significant bias when 
only high-mass subhalos are considered. In an earlier 
work Diaferio et al. (1999) found a bias for blue galax- 
ies of about 1.5 — 2, but no bias for red galaxies. Their 
interpretation was that the bias is due to the fact that 
blue galaxies are not in equilibrium. The velocity bias 
can be explained by the fact that slow subhalos are much 
less common, due to tidal disruption early in the merging 
process, before the cluster was virialized. Halos with low 
relative velocities can merge shortly after entering the 
cluster, thus decreasing the number of small subhalos 
with low velocities; after virialization, mergers are sup- 
pressed (Gnedin 2003). Large subhalos, > 10 -3 M c i uster , 
are more likely to merge with the central halo, rather 
than in subhalo-subhalo mergers (Angulo et al. 2008). 
Thus, in the central part of the cluster there is also a 
decrease in the number of massive subhalos with low ve- 
locities. This may explain why the value of b increases, 
going from the outer region towards the cluster center. 
The baryon content of halos may also affect the velocity 
bias: Due to ram-pressure stripping of galactic gas, which 
is more effective in the central cluster region (e.g., Arieli, 
Rephaeli, & Norman 2008), the subhalo mass is more 
easily reduced by tidal disruption. The deduced subhalo 
mass function is reduced relative to that deduced from a 
corresponding DM-only simulation (Dolag et al. 2008). 

The current data only allow determining an overall typ- 
ical value of /?dm; its best-fit value is in general agreement 
with the /3 ga i deduced from the data, and the value of 
/3dm found in simulations. As we noted, our assumption 
of eq. [3] is likely valid at large radii, which implies that 
the procedure specified in §[2]is more reliable there. Host 
et al. (2009) used the X-ray temperature as a surrogate 
measurement for deriving the DM velocity anisotropy 
profile at intermediate radii by assuming a connection 
between T gas and Tdm- They applied their analysis to 
11 low redshift and 5 intermidiate redshift clusters, in- 
cluding A1689. Their deduced values for /3dm are in very 
good agreement with our best-fit value. Together the two 
methods enable estimating the DM velocity anisotropy at 
all radii. Our approach here is less prone to the substan- 
tial uncertainty inherent in an attempt to determine the 
DM velocity anisotropy from a similar treatment which 
is based on IC gas properties. A possible difficulty with 
the latter approach may be due to a nonthermal pressure 
component that could appreciably affect results inferred 
from a solution to the hydrostatic equilibrium equation 
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(e.g., Molnar et al. 2010). 

It can be seen from figure |3] and @] that at r ~ 700 h -1 
kpc ~ ^r v i r there is a change in the relation between the 
dynamical properties of the DM and galaxies. Carlberg 
(1994) showed analytically that if one assumes a power- 
law mass density profile and isotropic velocity dispersion, 
this yields a power-law profile for the velocity dispersion. 
He also found that in general a cooler tracer has a steeper 
density profile and is more centrally concentrated. We 
find that at small radii, r < |r V j r , / co n < 1, implying 
that DM is cooler than galaxies. Indeed, at these radii 
the DM profile is steeper (see figure [4} and more concen- 
trated (see figure [3]) than that of galaxies. The subhalo 
distribution would naively be expected to be closely re- 
lated to the distribution of galaxies. In fact, galaxies and 
subhalos represent different populations and are not di- 
rectly comparable since subhalo masses are more strongly 
affected by tidal stripping than galactic baryonic matter. 
It was shown in simulations that subhalos within 0.3r v ; r 
typically lose more than 70% of their mass during the 
merging phase, while subhalos at r > 0.5r v i r typically 
lose only < 40% of their mass (Nagai & Kravtsov 2005, 
hereafter NK05). 

It was previously inferred from simulations that there 
is a spatially anti-biased subhalo distribution, in the 
sense that the DM subhalo profile has a larger core ra- 
dius than the background DM profile (DMS04; Gao et al. 
2004a) . DMS04 interpreted their result as due to missing 
half of the halo population caused by the known numer- 
ical overmerging problem. However, NK05 used a high- 
resolution cosmological cluster simulation to show that 
the subhalo radial distribution is significantly less con- 
centrated than that of DM due to tidal stripping, rather 
than this being a numerical artifact. They demonstrated 
that the radial bias disappears almost entirely if subhalos 
are selected using their mass or circular velocity during 
the merger phase. A similar result was obtained by Gao 
et al. (2004b) who found that defining the subhalo pop- 
ulation by requiring a minimum circular velocity gives 
a subhalo distribution which is more concentrated than 
when selection is based on a minimum mass. This re- 
sults from the fact that the subhalo distribution is easier 
to track by using the maximum circular velocity values, 
since mass loss is also accompanied by a decrease in the 
maximum circular velocity, but the decrease is slower 
than in the mass (NK05). 

In many previous studies the observed galaxy distri- 
bution was taken to be cuspy and similar to that of DM 
(e.g., DMS04; Gao et al. 2004a; NK05), using the galaxy 
distribution from Carlberg, Yee, & Ellingson (1997, here- 
after CYE97), who co- added observations of 14 clusters 
at various redshifts. The superposed sample contained 
1150 galaxies, including background galaxies. In the L09 
analysis of A1689, 500 cluster members were identified 
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from spectroscopic data, and about 1900 cluster mem- 
bers from photometric data. L09 showed that the galaxy 
density profile is best fitted with a core profile; although 
a cuspy profile was found to be acceptable, the deduced 
values of the scale radius and the power-law index were 
questionable. Indeed, the cuspy profile in CYE97 gave 
a better fit with a higher power-law index (preferring 
Hernquist (1990) rather than an NFW profile), yielding a 
quite high value of the scale radius, r s — (0.66±0.09)r 2u o- 
Moreover, Adami et al. (1998) examined a sample of 62 
clusters and found that most galaxy density profiles are 
better fitted with a cored rather than a cuspy profile, 
though for individual clusters the preference for a cored 
profile is rarely significant at the 90% confidence level. 
When Adami et al. (1998) composed a superposed sam- 
ple they obtained a clear preference for a core (King) pro- 
file (at more than 95 % confidence). Adami et al. (1998) 
showed that CYE97, as well as Beers & Tonry (1986), 
obtained a cuspy profile due to selection bias caused by 
not taking into account the effect of elongation. In some 
recent simulations (e.g., Saro et al. 2006) that included 
gas cooling and star formation, a detailed treatment of 
stellar evolution and chemical enrichment, as well as SN 
energy feedback in the form of galactic winds, the galaxy 
density profile has a shallower core, quite different from 
the cuspy DM profile at small radii, r < 0.4r 2 oo- 

Finally, we note that our approach of inferring /? ga i 
using the Jeans equation, which is derived from the col- 
lisionless Boltzmann equation, is not fully self-consistent 
since at small radii, r < 4r V i n we find possible deviation 
from a fully collisionless behavior. In addition, for sim- 
plicity we assumed that the cluster is spherically sym- 
metric, though there are claims that this cluster has a 
significantly triaxial shape (Oguri et al. 2005; Morandi 
et al. 2010). However, biases due to the triaxial nature 
of the cluster should affect the DM and galaxies fairly 
similarly since they have a similar distribution. More 
importantly, it should be kept in mind that the results 
reported here are based on a comprehensive analysis of 
only one cluster. Obviously, the results should be viewed 
as preliminary until reproduced by a similar analysis of 
a sufficiently large sample of clusters. We plan to do so 
with a larger sample of relaxed X-ray selected clusters in 
the CLASH program^. 
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